Transport through a molecular quantum dot in the polaron crossover regime 
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We consider resonant transport through a molecular quantum dot coupled to a local vibration 
mode. Applying the non-equilibrium Green function technique in the polaron representation, we 
develop a non-perturbative scheme to calculate the electron spectral function of the molecule in the 
regime of intermediate electron-phonon coupling. With increasing tunneling coupling to the leads, 
correlations between polaron clouds become more important at relatively high temperature leading 
to a strong sharpening of the peak structure in the spectral function. The detection of such features 
in the current-voltage characteristics is briefly discussed. 

PACS numbers: 73.23.-b, 71.38.-k, 73.63.-b, 85.65.+h 



The possibility of using the degrees of freedom of 
molecules contacted to metallic leads for nanoelectron- 
ics has been pointed out in the last decade [H, 0, 01, 0] . 
Molecular devices typically consist of single molecules, 
carbon nanotubes,..., connected to leads, where vibra- 
tional degrees of freedom and electronic interactions op- 
erate. Their characterization via transport remains an 
experimental and a theoretical challenge. Depending on 
the strength of the electron-phonon coupling, various ap- 
proximation schemes have been developed for weak cou- 
pling 0, 0, 0] and for the intermediate to strong cou- 
pling, where the electrons of the molecular quantum dot 



are dressed by a polaron cloud 

H S El EI El E3, 

EIE1E1E3. The purpose of the present paper is to 

go beyond the existing polaron approaches to molecular 
transport, allowing a better description of the intermedi- 
ate electron-phonon coupling regime. While conventional 
perturbation theory is typically formulated with the help 
of (single) phonon Green functions, for the intermediate 
coupling regime a more natural choice is to introduce a 
Green function which describes collective phonon excita- 
tions associated with the polaron cloud. We show that 
the electron spectral function of the molecular quantum 
dot can be evaluated systematically by diagram dressing 
of the polaron Green functions. In the polaron crossover 
regime, where the coupling to the leads becomes compa- 
rable to the phonon frequency and/or the strength of the 
electron-phonon coupling, the phonon sidebands in the 
spectral function are found to remain sharp even at high 
temperature. This implies that at intermediate coupling, 
phonon features for the current through the molecule 
could be more pronounced than expected. Pauli block- 
ing of linewidths in the strong coupling regime and at 
low temperature has recently been discussed [T3 |. 



We consider a generic model which couples a single- 
electron level to a local vibration mode with frequency Q. 
The Hamiltonian of the molecule coupled via tunneling 
to two metallic leads is given by H = H m + H\ ea ds + Ht, 



with (h = 1): 



~e -g(b + b r )] dU + Sltfb, 



(1) 



where d (b) are annihilation operators for the dot elec- 
tron (phonon), £q is the bare dot level, and g is the 
electron-phonon coupling. The lead Hamiltonian is 
c ]k c jk, where c jk is the electron oper- 
ator for lead j = I, r. The tunneling Hamiltonian is 
Ht = J2jkTjk.Cj k d + h.c, where Tjk is the tunneling 
amplitude between lead j and the molecule. For sim- 
plicity, we ignore the spin degree of freedom assuming 
a resonant tunneling situation in the Coulomb blockade 
regime, far from the Kondo regime. The bias voltage 
across the junction, V — Vi — V r , is imposed by shifting 
the chemical potentials of the leads, fij = eVj. Through- 
out the paper, we consider the symmetric case, Tjk = 
assuming the electron-hole symmetry to be held in the 
presence of phonons (see below) . In the steady state, the 
current through the molecule can be expressed as [18j ]: 



I{V) 



w)-/ r (w)]r(w)yl(w), (2) 



where A(ui) is the Fourier transform of the elec- 
tron spectral function of the dot, A(t — t') — 
-i({d(t), dt(f )}>, f l/r (u) = f(u T eV/2) are Fermi dis- 
tributions in the leads, f(ui) = (e' 3 " + l) , and T(uj) — 
2nJ2k \1~k\ 2 S (oj — £fc) is the tunneling half-width. 

We eliminate the electron-phonon coupling term in 
H m by using the polaron unitary transformation U [l9| : 
H = U\HU = H m + H leads + H T with U = e ~ ia P d ' d , 
p = i (w — 6), and a = g/fl. The diagonalized Hamilto- 
nian of the molecule reads: H m = ed^d + tlb'b, where 
e = €q — g 2 /CI is the renormalized level energy. The as- 
sumed above electron-hole symmetry implies that e = 
which can be achieved by applying a gate voltage. The 
transformed tunnel Hamiltonian becomes dressed by the 



vibrations, Ht = J^ik^k 



D 



where D = dX 



and X 



-lap 



is the polaron cloud operator. This 



form of Ht suggests an effective phonon-mediated cou- 
pling between the dot and lead electrons. In the polaron 
representation, the spectral function reads A(t — t') — 
-i({D(t), Dt(t')}). 

We thus introduce the Keldysh Green function for the 
polaron, Q ss ' (t - t') = -i (T K { D s (t)D^ s ' (£')}), where 
s, s' = {1, 2} stand for the forward (1) and backward (2) 
branches of the Keldysh contour, and Tk is the Keldysh 
time-ordering. Averaging over the leads results in: 

G ss \t - t') = -i (T K {S T D s (t)D^(t')}) , (3) 

where ( ...) = is a trace over the states of the molecule 
only, and the time-dependence of D(t) is given in the 
interaction picture with respect to H m . In Eq. (|3|), the 
effective S-matrix operator reads: 



St = Tk exp < — i 



+oo 



dTdT'W(T,T') \ , (4) 



W(r,r') = J2 (r)Sr' (r - r')D s ' (/) , (5) 

ss' 

where the electron tunneling self-energy Eo is built up 
from the Keldysh Green functions for the uncoupled 
leads, <$'(r), S^'(r) = (-1)*+°' £. fe |T fe | 2 gg(r), and 
can be expressed in terms of T(oj) and fj(uj) by virtue of 
E<'>(uO = ±iT( W )£j /i(±o;). 

We now proceed with the perturbation theory in the 
lead self-energy. The main challenge resides in the fact 
that because of the multi-phonon operator factor X ac- 
companying the d operator in St, the Wick theorem 
is not applicable to express Q in terms of the electron 
and phonon self-energies. It is convenient to introduce a 
Keldysh Green function describing the polaron cloud: 



A s 



(*)}, 
(6) 

with $ ss '(t) = (1/2) (Tjf (p*{t) -p s '(0)j ) • A typical 

approximation [ll|, [l2|, [3, [2(| , which is called the "single 
particle approximation" (SPA) in Ref. LLJ, consists in 
factorizing the electron Green function as 



g ss ' (t - if) ~ G ss ' (t - if) A ss ' (t - t') 



(J) 



where G obeys the Dyson equation (in the Fourier repre- 
sentation) G(w) = G(lu) + G{uS) E (w) G(w), with G the 
bare Green function of the isolated dot (in the absence 
of phonons). This decoupling of electron and phonon dy- 
namics can be justified in the anti-adiabatic regime [2l|, 
where f2 (and/or g) is large compared to the dot tunnel- 
ing rates. In this regime, the phonon mode is fast and 
reacts instantaneously (on the time scale of T _1 ) to the 
presence/absence of an electron on the molecule. Eq. ([7]) 
leads to the following expression for the spectral function: 



where A>(t) = A<(-i) = e^ 2 *^*) with 

$>(f) = (2JV + 1)(1 - cos fir) + i sin fir , (9) 
and N — (e^ n — l) . The Feynman diagrams associ- 
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FIG. 1: The Dyson series for Q corresponding to the 
single-particle approximation, Eq. 0. Solid lines 

represent the bare electron propagator G, dashed lines 
the self-energy Eo (lead electrons), and dotted lines 
indicate the phonon-cloud dressing propagator A. 

ated with the SPA are depicted in Fig. [2 note that in 
each order the phonon cloud connects only the end points 
of the full fermionic lines. In [3, [H| > the equation of mo- 
tion method was used to go beyond this approximation. 
Below, we compare the SPA with our own scheme, la- 
beled as the "nearest-neighbor crossing approximation" 
(NNCA), which incorporates the effects of the polaron 
cloud in the vertex part of the Green function. 

First, to zeroth order in the tunneling rates, the 
Keldysh Green function for the polaron is given by an 
expression similar to Eq. (|7|), 



G 6 s (t - if) = G ss (t - O A ss ' (t - t'), 



(10) 



A SPA {t) = G > (t)A > (t)-G < (t)A<(t), 



(8) 



but with the bare propagator G. In this so-called atomic 
limit, Qo describes the localized electron state dressed by 
vibrations. This term corresponds to the first diagram 
in the series shown in Fig. [TJ Next, to first order in 
the tunneling rate, one obtains: Qf s (t — if) = Qf s (t — 
t')A ss '(t-t'), with 

QT\t -t')= f d Tl dr 2 GP(t - n) [A SS2 (t - T2T 1 

x E SlS2 (n - r 2 ) [A SlS '(n - O] ^ G s 2S '(r 2 - t') , (11) 

where we have introduced the dressed self-energy 
E SiS2 (t) = E^ lS2 (r) A SlS2 (r). We note that in Eq. (TiT]). 
the phonon-cloud propagator A^ 1 plays the role of a 
"phonon undressing" factor. 

Although Eq. (fTT|) looks rather complicated, it can be 
drastically simplified if one assumes that Re A 3> ImA, 
which is equivalent to saying that A K S> A R — A A , where 
A R I A I K is the retarded/advanced/Keldysh component of 
A. Neglecting ImA turns out to be a good approximation 
for the case of intermediate electron-phonon coupling, 
a < 1. This is illustrated in Fig. where it is shown that 
the real part of A dominates with decreasing electron- 
phonon coupling a. As seen from Eq. ([5]), this result 
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FIG. 2: Ratio between the imaginary and real parts of 
A> as a function of time t for a — 0.8 (dotted line), 0.6 
(dashed dotted line), 0.4 (dashed line) and 0.2 (full 
line). 

does not depend on temperature. Noticing that 

A ss ' = (1/2)A K + (1/2) (A ss ' - A 55 ') , (12) 

with s — 2(1) for s — 1(2), and neglecting the last term 
in the above expression which is proportional to ImA, 
one finds from Eq. pip for the retarded component: 

Q?\t -t') = J d Tl dT 2 g R (t - n) A- 1 (t - t 2 ) 

x e«(ti - r 2 ) A-^n - 1>) g R {r 2 - 1>) , (13) 

where for notation convenience we have introduced 
A(t)=A*(t)/2. ^ 

Now, going to higher order terms (n > 1) of per- 
turbation theory in the lead coupling, we replicate the 
structure of the first order diagram [n = 1) keeping only 
the nearest-neighbor crossing lines (A" 1 ) connecting the 
dressed electron propagators (Go) and the dressed self- 
energies (E). By doing this, we take into account cor- 
relations between phonon clouds corresponding to two 
sequential tunneling of the lead electrons, which is hap- 
pening on the time scale of the order of r . Noticing 
that the general solution for the Green function can be 
written as 

g R (t - o = g R (t - o + a(* - o g R (t - f) , (w) 

the n-order term for g R in the NNCA reads (n > 2): 

G?(t - = / dndTi ... dr 2n - n ) 

x \-\t - r 2 ) E fl (n - r 2 ) A- X (n - r 3 ) ... 
... E fl (r 2n _! - r 2 „) A _1 (t 2 „_i - t') g o fl (r 2n - . (15) 

The Feynman diagrams associated with this approxima- 
tion scheme are illustrated in Fig. [3l Going to the Fourier 
representation and summing up all the terms of pertur- 
bation theory, one finds: 

g R (uj) = g^{uj) u r (uj) g*{uj) , (ie) 



FIG. 3: Diagrammatic representation of the NNCA. As 
in Fig. [Tl dotted lines above the fermionic line represent 
the phonon-cloud dressing propagator A, while dotted 
lines crossing below the fermionic line represent the 
undressing corrections described by A -1 . 

where g R (oj) is the Fourier transform of 

g H (T) = \- 1 (T)g H (r), a?) 

satisfies the Dyson equation, 

n R (Lu) = T, R {u;) + Z r (uj) K r (lu) IFV) , (18) 

and K (u>) is the Fourier transform of 

K R (r)=X- 2 (r)g R (r). (19) 

Eqs. (fl4" |) -([T9" ]) constitute the basis of the NNCA approx- 
imation developed here. 

Although, after dressing, the retarded self-energy 
E«(t) = e(r)[A>(r)E>(T)-A < (r)E<(r)] becomes 
voltage-dependent, inspection of this expression in the 
wide-band limit, where the energy dependence of T(lj) 
is neglected, shows that for the intermediate coupling 
a < 0.5 the magnitude of T, R is only slightly decreased 
compared to the bare one, Eq 2 = — iT (see inset of Fig. 
[4}. As a result, in the wide-band limit considered here, 
the voltage dependence of the spectral function can be 
discarded. 

The difference between the NNCA and the SPA is hid- 
den in the kernel K R entering the Dyson equation (|18p . 
To show this we first notice that the expression for Q R , 
^(r) = A- 1 (r)e(r)[G>(r)A>(r)-G<(r)A<(r)], can 
be approximated by the bare Green function, g R w G R , 
which is valid up to the small term of the order of 
ImA/ReA. As a result, the kernel K R , Eq. HI]), can 
be viewed as the "anti-dressed" propagator, K r (t) w 
A _1 (r)G fl (r). This expression is different by the fac- 
tor A -1 from the corresponding expression in the SPA, 
K R PA — G R . This can be checked by replacing K R in 
Eq. (fT5)) by G R , which results in the SPA expression for 
the retarded Green function: 

g R PA (t-t')*\(t-t')G R (t-t'), (20) 

with G R (u>) = (u - e - E fl (w)) _1 . 

In Figs. [4] and [5] we plot the results for the spec- 
tral function A(u>) = — 2lmg R (uj) calculated within the 
NNCA by solving numerically Eqs. (T4"|) - (fll))) . For com- 
parison we also show the results corresponding to the 



4 



SPA based on Eq. ([8]). In all cases we take a = 0.4 
and e = 0. In the low temperature regime (Fig. [4]), one 
sees that the NNCA predicts a slightly sharper central 
peak and more pronounced phonon sidebands than ex- 
pected in the SPA at relatively high transparencies. In 
this regime, the peak sharpening can be attributed to the 
suppression of the tunneling self-energy by the phonon 
cloud, which pins the dressed electron on the dot (inset 
of Fig. |4|). Because of the low temperature, only the 
single-phonon peaks are clearly visible. With increas- 




FIG. 4: The spectral function at T/Q = 0.4 calculated 
within the NNCA (solid line) and SPA (dashed line) for 

a = 0.4, 2T/Q = 0.5, e = 0, and V = 0. Inset: 
Renormalized lead self-energy for (from top to bottom) 
a = 0.2, 0.4, and 0.6. The other parameters are the 
same. 

ing temperature, the occupation probability of phonon 
sidebands increases, and more satellite peaks appear in 
the spectral function (Fig. [5]). Although for small T, as 
expected, the NNCA and SPA results are hardly distin- 
guishable (see inset of Fig. [5]), the difference is clearly 
seen for the higher transparency (the main plot of Fig. 
[5]). According to the NNCA, all excited satellite peaks re- 
main sharp and are rather robust with respect to increas- 
ing r. This is the regime where the SPA breaks down be- 
cause the tunneling rate (2r) becomes comparable to Q, 
In this transition region, increasing the temperature acts 
in favor of the electron-phonon interaction which tends 
to localize the dot electron, leading to a well pronounced 
peak structure in the spectral function. With increas- 
ing the T, such that 2r > Q, we expect that the NNCA 
will not be sufficient, and the higher order correlations 
between the polaron clouds should also be included. 

Finally, we remark that the effect of suppression of the 
tunnel broadening in the spectral function at high tem- 
perature can be observed directly in the measurement of 
the differential conductance keeping the leads at very low 
temperature (T <C Q), so that dl/dV oc A(eV), but vary- 
ing the local temperature of the molecule, for instance, 
by heating it with a laser. 

We thank D. Feinberg for useful discussions. 
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FIG. 5: Same as Fig. g]but for T/Q = 4 and 
2r/0 = 0.8. Inset: same as the main figure but for 
2T/Q = 0.2. 
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